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ABSTRACT 

We demonstrate a novel mechanism for the formation of topological defects in a first 
order phase transition for theories in the presence of small explicit symmetry breaking 
terms. We carry out numerical simulations of collisions of two bubbles in 2+1 dimensions 
for a field theory where U(l) global symmetry is spontaneously as well as explicitly broken. 
In the coalesced region of bubble walls, field oscillations result in the decay of the coalesced 
portion in a large number of defects (e.g. ten vortices and anti- vortices). We discuss the 
implications of our results for axionic strings in the early Universe, for baryon formation 
in quark-gluon plasma, and for strings in liquid crystals in the presence of external electric 
or magnetic field. 



Production and subsequent evolution of topological defects has been of considerable 
interest for particle physicists in the context of the early Universe |]T|]. Similar tech- 
niques have also been used to study baryon formation during hadronization of quark-gluon 
plasma (QGP) in heavy ion collisions H- Study of topological defects has, of course, 
been possible in most detailed way only in condensed matter systems where they can be 
experimentally studied @]. 

The aforementioned defects correspond to a spontaneous breakdown of a symmetry. 
However, there are many situations when the symmetry is also explicitly broken. In 
particle physics, the Pecci-Quinn scheme for solving the strong-CP problem of quantum 
chromo dynamics leads to the presence of an explicit symmetry breaking term and, con- 
sequently, to axionic strings 0. The Skyrmion picture of baryons in the context of chiral 
models is another example where explicit symmetry breaking terms are needed to incor- 
porate a non-zero pion mass 0. In condensed matter, liquid crystals provide a simple 
example of such systems where the presence of external electric or magnetic fields induces 
explicit symmetry breaking terms 

The study of formation of topological defects in such systems is therefore important as 
it has implications for a diverse set of phenomena. It has recently been argued that explicit 
symmetry breaking can lead to a four fold enhancement in the production of baryons in 
QGP 0. These arguments were largely qualitative and did not depend sensitively on 
the order of the phase transition. It was argued in that a similar enhancement should 
occur for other topological defects as well. 

In this letter we demonstrate a new mechanism for the production of topological defects 
for systems with explicit symmetry breaking and with a first order phase transition, where 
phase transition proceeds via bubble nucleation. This mechanism leads to a much stronger 
enhancement in defect production, and results from a combination of the effects discussed 
in as well as effects coming from the large field oscillations in the region of coalesced 
bubble walls. The net result is that wall oscillations decay by producing a large number 
of vortices and antivortices. For example in one simulation, we found 5 vortices and 5 
antivortices produced in a single two-bubble collision. 

We adopt the same numerical technique as used in previous simulations of vortex 
formation via bubble collision, see [^. We will study vortex formation in 2-1-1 dimensions 
in a field theory system described by the following Lagrangian 

where (f) and 6 are the magnitude and the phase of the complex scalar field $ ($ = (f)e^^). 
This Lagrangian describes a theory where the U(l) global symmetry is spontaneously 
broken, except for the presence of the last term which breaks this U(l) explicitly. When 
K is zero, this theory allows for the existence of a cylindrically symmetric vortex which 
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is a solution of the time independent field equations. For a non-zero k, the vortex looses 
azimuthal symmetry and is not a solution of the time independent equations of motion 
any more. 

For K = 0, the process of vortex creation via bubble nucleation has been described in 
0. At zero temperature, bubbles of true vacuum nucleate via quantum tunneling in the 
background of the metastable vacuum with = 0. These are described by the bounce 
solution which is an 0(3)-symmetric, least action, solution of the Euclidean field equation 



^ + - vw = (2) 

dr^ r ar 

where V{(j)) is the effective potential in Eq.(l) (with k = 0) and r is the radial coordinate 
in the Euclidean space. In the Minkowski space the profile of the nucleated bubble is 
obtained from the solution of Eq.(2) by putting t = 0. This bubble then evolves according 
to the classical field equations obtained from the Lagrangian in Eq.(l) in Minkowski space, 

°^^ = -^>^ = l,2 (3) 

where $ = $i +i$2, ^ is the d'Alembertian and time derivatives of fields are set equal to 
zero at t = 0. In a phase transition, 6 varies randomly from one bubble to another. These 
bubbles expand and vortices form at the junction of three or more bubbles if the phase 
6 traces a nontrivial winding in that region. This is the conventional Kibble mechanism 
of defect formation |^ which leads to the probability of vortex formation for 2 space 
dimensions equal to 1/4 per bubble @]. 

We wish to study the case when k is non-zero. First we briefly recall the physical 
picture described in 0. Consider a two bubble collision with the phase 9 in the two 
bubbles taking values vr + a and ir — a where a is small. As the bubbles collide, 6 in the 
coalesced portion will assume value tt due to the geodesic rule (essentially to minimize 
energy) and will keep evolving towards zero inside the bubbles. It is then easy to see that 
this leads to a winding one getting created on one end of the coalesced wall and winding 
minus one on the other end [^. It was argued in that this leads to roughly four fold 
enhancement in the number density of vortex production per bubble. 

However, it turns out that the actual dynamics of vortex creation has much richer 
structure, especially for a first order phase transition. As 6 in both the bubbles evolves 
towards zero, the coalesced portion of the walls undergoes large oscillations. Such oscil- 
lations have been described in ^ where it was shown (for the case of subcritical bubbles) 
that, as undergoes large oscillations, it passes through = forcing 9 to change to 

e + TT. 

This flip in the orientation of $ has very important effects on the process of vortex 
formation. The evolution of 9 towards zero inside the bubbles tends to create a winding 
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one near one end and anti winding near the other end inside the coalesced region. The 
flip in the orientation of $ in the central region completes these windings and results in 
the nucleation of a vortex-antivortex pair in the coalesced wall. 

This explains the formation of the first vortex-antivortex pair. Subsequent pairs are 
created due to the oscillation of 9 about = 0. As ^ in the two bubbles evolves towards 
zero, it overshoots the true vacuum (i.e. = 0). It is easy to convince oneself that this 
evolution of 9 in the two bubbles, combined with the flipping of the orientation of $ due 
to wall oscillations, will result in the creation of another vortex-antivortex pair. This 
process continues as long as 9 and oscillate about 6* = and = respectively. 

It is important to realize that if the oscillations do not have large enough amplitudes, 
and associated with it appropriate 9 oscillations, then vortices are not nucleated. We see 
this in our simulations where many oscillations of the wall, and that of 9, may pass by 
before a given pair gets nucleated. Also, the sequence of vortices and anti-vortices created 
on one end of the coalesced region is quite arbitrary, depending on the details of and 
9 oscillations; though, the net winding number is always zero. This implies that the 
annihilation of vortex-antivortex pairs may be very ineffective. 

We now proceed to describe our numerical results. We find the bubble profile by 
solving Eq.(2) for k = in V{(f)). This bubble profile will be an approximate solution 
of the field equations obtained from Eq.(l) for small non-zero values of k, and provides 
an adequate starting point as bubbles collide only after undergoing large expansions. We 
choose the values of the parameters as, e = 0.2, 0o = 4.0, A = 4.0. We use natural system 
of units with fi = c = 1 and measure lengths in terms of rn]^ ^ where mu is the Higgs 
mass in the broken phase for the case when k is zero. = 0.12 for the above set of 

parameters. We have studied a range of values of n. Large values of k,{> 0.12) do not 
give any vortex formation as 9 in bubbles rolls down and settles to zero before bubbles 
can effectively coalesce. For all other values of k vortices form (with smaller k leading 
to vortex formation at a later stage). The figures shown in this paper correspond to the 
choice of K = 0.06. 

Following the techniques developed in [01 , we study the case of a two bubble collision 
by prescribing the nucleation centers for the two bubbles. This amounts to replacing a 
portion of the false vacuum region (with $ = 0) by the profiles of two bubbles. The 
field configuration is then evolved by using a discretized version of Eq.(3). Simulation is 
implemented by using a stabilized leapfrog algorithm of second order accuracy in both 
space and time. We use a 1000 x 1400 lattice with the lattice spacing in spatial directions. 
Ax, equal to 0.104 (in the units of mjj^) and lattice spacing in the temporal direction. 
At, equal to Ax/^/2. With these values, the evolution was completely stable and energy 
was conserved within a few percent during the simulation. For details of the numerical 
technique, see 0. Simulations were carried out on a HP- 735 workstation at the Institute 
of Physics, Bhubaneswar. 
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Bubble centers were chosen to lie at the y boundaries of the lattice (and at the midpoint 
of the X axis) so that the initial bubble profiles are that of half bubbles. We use free 
boundary conditions. Fig. la shows the plot of $ and the contour plot of (j) for the initial 
field configuration of one of the bubbles. Values of 9 for the two bubbles are taken as 
0.56 and 0.44 radians (for the lower and the upper bubbles respectively) and are uniform 
inside each bubble. This leads to development of a region of 6* = vr in the region where 
bubbles coalesce. Fig. lb shows the plot of $ at an intermediate stage. The bubbles 
have significantly coalesced and 6 inside the bubbles has started rotating towards zero. 
The rotation of 9 is smaller near the bubble walls due to the dependence of the explicit 
symmetry breaking term in Eq.(l) on (f?. [This is the reason that, even if the initial values 
of 6 in the two bubbles are not close to vr, a region of 6* = vr developes in the coalesced 
region and vortices are produced.] 

Fig.lc and Fig. Id show the configurations after the first pair has been nucleated. The 
plot of $ clearly shows that 9 has overshot the true vacuum (^ = 0) . Afterwards, 9 starts 
climbing towards tt first and then again rolls back towards zero. As described earlier, this 
will cause creation of subsequent pairs for appropriate (f) oscillations. Fig.le and Fig. If 
show the plots at a stage when there is a total of 10 vortices and antivortices. Out of 
these, there are two groups, containing three vortices each (as confirmed by detailed plots 
of $ of these regions) which are not well separated. Over all there are at least 6 vortices 
and antivortices which are well separated. Note that, due to the presence of explicit 
symmetry breaking term, the profiles of these vortices are highly deformed as shown by 
the contour plots. 

In conclusion, we have demonstrated a new mechanism for the formation of topological 
defects in the presence of explicit symmetry breaking which may dominate over other 
mechanisms of defect production. A somewhat modified version of this mechanism (due 
to the absence of coalesced portion of bubble walls) may also be applicable for the case 
of second order phase transitions. [In this context we mention that in [0 it was found 
that the number of vortices produced was roughly twice of the estimate based on the 
Kibble mechanism (though many pairs annihilated quickly). In view of our results in this 
paper, the excess production in [^ suggests that a nontrivial dynamics of 9 coupled with 
oscillations may contribute to vortex production even when explicit symmetry breaking 
is absent.] The most interesting aspect of this mechanism is that it is literally a pair 
creation process, though still governed by classical equations of motion. In this sense, it 
resembles the pair creation of vortices in the flow of superfiuid ^He through a small orifice, 
as discussed in [ffUl, though actual mechanisms are completely different. In a subsequent 
paper we will present the study of full phase transition by nucleating large number of 
bubbles [^ . 



Implications of these results are many. Using the ideas described above, it is possible 
to argue that, in two bubble collisions in 3+1 dimensions, field oscillations should lead to 
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string loops being emitted out from the coalesced region. For axionic strings in the early 
Universe, earlier studies have assumed that the formation mechanism is the same as for 
other cosmic strings, namely, via the Kibble mechanism [^]. Above discussion shows that 
the dominant mechanism may be via the mechanism discussed in this paper, at least for a 
first order phase transition. Therefore, the final distribution of axionic strings, and hence 
the frequency distribution of emitted axions, can be drastically different from what is 
conventionally taken. Also, as now one expects small string loops to be produced, axionic 
domain walls may not survive for long. This may make a larger class of axionic models 
viable. 

For the case of liquid crystals in the presence of electric field, our results suggest 
that, instead of long strings, small string loops should form in the coalesced region of two 
bubbles. However, in this case, the dynamics of string formation via this mechanism may 
be completely dominated by the presence of damping terms. The mechanism discussed in 
this paper should also be applicable to the production of other defects, though the details 
of the mechanism will depend on the type of defect (and dimensionality of physical space) 
under consideration. Especially important is the production of baryons (in the Skyrmion 
picture) in quark-gluon plasma . It is clearly important to investigate the enhancement 
(which may be very large, as suggested by our results) expected in baryon production due 
to explicit symmetry breaking terms if the phase transition is of first order. 
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FIGURE CAPTIONS 



Fig.l: (a) Plot of $ and contour plot of (f) for the initial configuration of one of the 
bubbles with center at the boundary. For all $ plots, the orientation of the arrows from 
positive X axis gives the phase of $ while the length of arrows is proportional to 0. (b) 
Plot of $ for the coalesced region at t = 55.2 showing that 9 has significantly rotated 
towards zero in bubble interiors, (c) Contour plot of at t = 62.6 showing a vortex- 
antivortex pair, (d) Winding numbers of the vortex and the antivortex are clear from the 
$ plot, (e) Contour plot of at t = 84.7 showing ten vortices and antivortices. There are 
two groups of three overlapping vortices each, one near x = 34 and the other near x = 
132. These groups have net windings of +1 and —1 respectively, (f) $ plot for a portion 
of lattice showing winding numbers of at least 2 vortices and 1 antivortex which are well 
separated. The winding +1 region near x = 34 actually consists of closeby configurations 
of two vortices and one antivortex. 
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